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We investigate low-energy excitations of the one-dimensional half-filled S'C/(4) Hubbard model 
with an attractive on-site interaction (7 < using the density matrix renormalization group method 
cis well as a perturbation theory. We find that the ground state is a charge density wave state with a 
long range order. The ground state is completely incompressible since all the excitations are gapful. 
The charge gap which is the same as the four-particle excitation gap is a non-monotonic function of 
U , while the spin gap and others increase with increasing |i7| and have linear asymptotic behaviors. 
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PACS numbers: 71.10.Fd, 71.10.Pm, 71.30.+h, 05.10.Cc 

The Hubbard model is one of the most classic mod- 
els for strongly correlated electronic systems and has at- 
tracted long-term interest since the pioneering work in 
1960s[i|. By taking the on-site Coulomb interaction into 
account, it well explains the puzzle that some materials 
with half-filled band are insulators. However, this model 
which consider only a single band, on-site Coulomb inter- 
action and the nearest-neighbor hopping is often thought 
being oversimplified. To account for other features be- 
yond the Mott physics Q, there are two kinds of natu- 
ral extensions: one is to incorporate the orbital degree 
of freedom which may be called multi-band Hubbard 
model [3!], and the other is to consider the hopping and/or 
the Coulomb interaction with longer ranges. 

In the last several years, ultra-cold atomic experi- 
ments evoke systematic studies of correlation effects in 
the optical lattice systems where interactions are tun- 
able through Feshbach resonance. The Hubbard model 
becomes again an appropriate one to envisage some rele- 
vant issues with both positive and negative interactions. 
Recently, fermionic atoms with higher spins are success- 
fully trapped into optical lattices [j] . This calls for a gen- 
eralization of the SU{2) Hubbard model into the SU{N) 
case t5l,l6|. In this paper, we study one-dimensional SU (4) 
Hubbard model which is represented as 
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U 
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where < > is a hopping matrix element, L the number 
of the lattice sites, a and a the spin indices taking — |, 



i and 2. 



Ci^icia-) denotes the creation(annihilation) 
operator of a particle with spin a at the site i, riia- — 
c\fjCia is the corresponding number operator and U is 
the on-site interaction. 

This model is not exactly solvable even in one di- 
mension in contrast to the SU{2) case[3, S]- Never- 



theless, some aspects of physical properties can be re- 
liably explored by some analytical approaches as well as 
numerical methods such as density matrix renormaliza- 
tion group (DMRG)[i, and Quantum Monte Carlo 
simulations [ill [l^ . Most recently, an 50(8) symmetry 
regime was proposed between < ?7 < 3i at half-filling 
by Assaraf et al[l^ with using a nonperturbative renor- 
malization group method and quantum Monte Carlo sim- 
ulation. They found that the low-energy spectrums are 
gapful in this regime. The similar results were also shown 
later by Szirmai and Solyom for the other N > 2 caseilS]. 
Those studies were concentrated on the repulsive case 
U > 0, but for the U < case few results are obtained 
so far. On the other hand, it is well known that the one- 
dimcnsional attractive half-filled SU(2) Hubbard model 
is described by a Luther- Emery liquid [lj|, in which the 
charge excitation is gapless, whereas the spin excitation 
is gapful. By the hidden SU{2) transformation the SU (2) 
Hubbard model with U can be mapped to the one with 
—U, while for the SU{4) case such a mapping does not 
exist so that one cannot obtain any insights into the low- 
energy properties through the mapping. In this paper, we 
will show that the SU{4) Hubbard model at half-filling 
with the attractive interaction belongs to a different uni- 
versality class from the SU{2) one. 

Let us start with a perturbation theory for the strong 
coupling regime. For this purpose, we rewrite the Hamil- 
tonian ^ as TL — TLt -\- Hu, where the hopping term 
Ht = ^tj^iai'^ia-'^i+icr + h.c.) IS regarded as a perturba- 
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tion and the on-site interaction Tiu = 
as the zeroth order Hamiltonian has highly degenerate 
ground states in which each site is either fully occupied 
by four particles forming a SU{4) singlet, or empty. Up 
to the second-order, the effective Hamiltonian is given by 
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where P is a projection operator which projects a state 
onto the subspace spanned by the ground states of 7i„, 
and rii = the number operator at the site i. The 

hopping term Ht hfts the degeneracy of W„ and gives rise 
to both the energy gain of ^ per site at half-fiUing as 
denoted by the first term of ^ and an effective repul- 
sive interaction between particles on the nearest-neighbor 
sites as denoted by the second term of The second 
term induces essentially a charge-density-wave (CDW) 
ground state with a true long range order such that ev- 
ery other site is fully occupied with an empty site in 
between, which is consistent with the mean field result 
for the weak coupling region @. Moreover, H^^j has two- 
fold degenerate ground states, each of which is a SU (4) 
singlet. 

An important question concerning the CDW ground 
state is whether it is metallic or insulating. To address 
this question, we need to examine charge excitations. 
The charge gap Ac is the energy difference of the lowest 
excitation in the spin singlet channel from the ground 
state as defined by = Ei{L,2L,Q) ~ Eq{L,2L,0), 
in which En{L, N, S) stands for the the n— th excita- 
tion energy in a spin-S* channel with L sites and N 
particles. And another interesting issue is to explore 
the relevance of the four-particle excitation gap A4 to 
Ac for the 5'C/(4) symmetry. A4 represents the en- 
ergy cost of adding four particles or holes into the sys- 
tems such that A4 = i [Eo{L,2L + A,0) + Eo{L,2L - 

4,0) -2Eo{L,2L,0)]. From H^^'^, one can easily find 

Ac = — IIj- based on the fact that the motion of four 
particles from one of fully occupied sites to its neighbor 
costs a minimal energy. Similarly, adding four particles 
to the system would gain the energy 6J7-|-2x (— while 

adding four holes costs — 6C/, thus one has A4 = — It 
turns out that the four-particle excitation gap is essen- 
tially the same as the charge gap and finite. 

However, the situation is completely different for the 
SU{2) case, where the ground state is metallic. To un- 
derstand this, one can write the effective Hamiltonian in 
the strong coupling limit for the SU{2) case as[l5|: 

i(7 

(2) 

This Hamiltonian distinguishes itself from Ti-lj^ with its 
last additional term which allows a "pair hopping" pro- 
cess between the two neighbor sites. Although the sec- 
ond repulsive term would be in favor of forming a CDW 
ground state with gapful charge excitations, the extra 
"pair hopping" term eventually destabilizes this CDW 
long range order resulting in gapless charge excitations. 
On the other hand for the 5J7(4) case, similar hopping 



term can only appear after calculating to higher order. 
Therefore, the perturbation theory in the strong coupling 
regime shed light on the different nature of the ground 
states between the SU{2) and SU{4:) cases. 

In order to go beyond the validity range of the above 
perturbation theory for a full exploration of the low- 
energy properties, we have performed systematic DMRG 
computations. There are actually some difhculties inher- 
ent to this model, such as a large number of degrees of 
freedom for each site and different edge states so that 
some measures are taken necessarily to reach sufficient 
accuracy in our computations. To calculate Ac, we have 
to use the periodic boundary condition (PBC) with the 
even number of sites because of multi-edge excitations 
in the SU{A) singlet subspace when open boundary con- 
dition (OBC) is imposed. For other gaps, we can effi- 
ciently expel the corresponding edge excitations by an 
OBC algorithm. In particular, one lattice site is added 
at each step and broken into two pseudo-sites. When 
the infinite system algorithm is conducted with the size 
of the supcrblock up to some odd number of sites L 
which are preselected, the sweeping procedure is per- 
formed for those L. The necessary extrapolations for 
the thermodynamic limit are finally made properly on 
the data with these preselected sites. In this case, we 
have to redefine the gaps correspondingly, for instance 
A4 = Eo{L, 2L + 6, 0) - Ea{L, 2L + 2,0)- 6C/, where the 
particle- hole symmetry is explicitly taken into account. 
In the strong coupling region, PBC is often used to iden- 
tify the bulk values of a gap rather than edge-excitation 
energy obtained under OBC. In our computations, t is 
set to be unit and 2000 states are kept for most cases 
and the maximal truncation error is the order of 10~^. 
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FIG. 1: Ground state energy correction per site from n^f^ 
(solid line) versus U in comparison with DMRG data (O)- 

Figure [1] displays the ground state energy correction 
per site arising from the hoping term for the thermody- 
namic limit. Both the numerical results and the pertur- 
bation theory show this energy increases monotonically 
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as \U\ increases and approaches to zero asymptotically. 
When — 2 < [/ < 0, the results of the perturbation theory 
deviates from DMRG ones, and the deviation becomes 
significant for small \U\. This is reasonable because in 
this region the hopping term Tit is no longer perturba- 
tivc. However, for U < —2, the perturbation theory pro- 
vides very good results which accurately agrees with the 
DMRG data. On the other hand, our DMRG calcula- 
tions with both OBC and PBC show that for a finite and 
even L, the ground state is unique and a iS'C/(4) singlet 
which belongs to the irreducible representation [1*^] [iq. 
In addition, slightly above the ground state there is one 
accompanied SU{A) singlet excited state, whose energy 
difference from the ground state diminishes as L — > oo. 
Therefore, one obtains two-fold degenerate ground states 
in the thermodynamic limit, which are consistent with 
our analysis based on Ti-eff These degenerate CDW 
ground states with the long range order result from the 
translational symmetry breaking. 

In FIG. m we show the DMRG results on the charge 
gap Ac and the four-particle excitation gap A4 for the 
entire range of ?7 < 0. First of all, one can see that 
both Ac and A4 are non-vanishing for all finite U < 
so that the ground state is insulating rather than metal- 
lic in contrast to the SU{2) case. Secondly, these two 
gaps behave non-monotonically with U. The maximum 
shown around U = —2 indicates a crossover region be- 
tween weak and strong interaction regimes. Third, the 
perturbation theory for the strong coupling regime pro- 
vides correctly the asymptotic behavior for large \U\ limit 
and shows a qualitative agreement with the DMRG re- 
sults in the strong coupling regime. The visible deviation 
from the DMRG results sets on at about U ^ —5 lower 
than that {U w —2) for the ground state energy correc- 
tion shown in Fig. 1. In the weak coupling regime, Ac 

and A4 shown by DMRG decrease with increasing U. Fi- 
(2) 

nally, while He/f can predict Ac — A4 only in the strong 
coupling limit, our DMRG calculations show that within 
the numerical accuracy Ac remains equal to A4 beyond 
the strong coupling regime. Although it is difficult from 
the DMRG calculations to obtain sufficiently accurate Ac 
for — 1 < J7 < yet, it is reasonable to conclude that A4 
is equal to Ac for all C/ < in the SU (4) case. 

Now we turn to the other three types excitations: the 
first one is the quasi-particle gap Ai for adding single 
particle or single hole to the system, the second one spin 
gap As corresponding to the excitation energy in the 
spin triplet channel from the ground state, and the last 
one two-particle gap A2 defined as energy cost when two 
particles or two holes are added to the system. While 
these three gaps together with Ac and A4 essentially in- 
volve all kinds of relevant excitations, they have signifi- 
cantly different behaviors. Since the ground state is the 
CDW state with the long range order, it is insightful to 
analyze those excitations in the Hartree-Fock (HF) ap- 
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FIG. 2: DMRG results for charge gap (+) and four-particle 
gap (O) versus U in comparison with perturbation ones (solid 
line). Error bars are smaller than size of symbols except for 
Ac at U = — 1 as estimated with keeping different states. 

proximation. The CDW state in HF approximation can 
be achieved by simply writing the on-site interaction as 
n-icrn^^' « {ni„)n^^> + nicr{n.^^') - {ni„){n^^') and assum- 
ing {riicy) = n -\- (— I)*(5n, where n is the average number 
of particles per site for each spin and 5n the correspond- 
ing order parameter. At half-filling, one has n — ^ and 
< (5n < i. By further introducing a;o- = C2ia- and 
bia — C2i+ia for cach sublattice of the bipartite lattice, 
respectively, and taking the Fourier transformation, then 
we can write down the Hartree-Fock Hamiltonian as: 

- -t^((f + e-^)a+>fe. + /i.c.) 

k(7 

+W5n 'Y^ia'l^aka - Ka^ka) + const. (4) 

k(7 

Diagonalizing this Hamiltonian, one can obtain two 
bands for each spin species a with the quasi-particle dis- 
persions = ±^A^ -f- 4^2 cos^ I where Ai — —3Udn. 
Moreover, one has Ac = A^ = 2Ai in the HF ap- 
proximation. These gaps can be then evaluated af- 
ter solving the following self-consistent equation 26n — 
(a^^aia) — {b'l^bic) for the order parameter 6n. In or- 
der to calculate the two-particle excitation gap, however, 
it is necessary to account for the particle-particle cor- 
relations into the HF approximation, which is nothing 
but the random phase approximation (RPA). For this 
purpose, one first constructs the basis with two-particle 
excitations from the ground state of Ti.^^ as follows: 

= «+ I*,) ^\{Ptm (5) 

ka 

where a^^ is an operator creating one quasi-particle with 
momentum p and spin a in the band and in \^ g) the 
band w~ is fully filled up by quasi-particles fi'^^ with the 
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momenta k and spin a. Then A2 can be obtained by 
diagonalizing H on the above basis ([5]). 
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FIG. 3: DMRG results for one particle gap Ai (x), two- 
particle gap A2(n), spin gap As (A), four-particle excitation 
gap A4 (O) and charge gap Ac(+) are shown as a function 
of U. The results of HF-RPA for Ai, A2 and A^ denoted by 
solid, dot-dashed, and dashed lines, respectively. Inset shows 
A/\U\ versus U for Ai (x), A2 (□) and A^ (A). 

As compared to the DMRG results, we found that HF- 
RPA can provide a qualitatively correct description for 
Ai, A2, and A^. Figure [3] shows the DMRG data on 
all five gaps as well as Ai,A2, and A^ from HF-RPA 
approximation for the region of — 2 < [/ < and the 
inset illustrates Ai, A2, and A^ by showing ratios for 
them over \U\ up to \U\ = 10. In contrast to A^ and 
A4 as seen from Fig. 2, Ai, A2, and Ag increase with 
increasing \U\ and become linear in large \U\ limit. It 
turns out that the relation Ac = A^, given by the HF 
approximation, is invalid for general U < 0. Moreover, 
it is unclear but beyond the present approaches whether 
there is a symmetry enlargement similar to the one pro- 
posed for the repulsive case[l3|. Nonetheless, HF-RPA 
presents precise asymptotic behaviors for Ai, A2 and 
Ag. In the weak coupling limit, the exponential open- 
ing of these gaps can be well reproduced from the solu- 
tion to the self-consistent equation 6n ~ and 
with taking into account the two-particle correlations. 
In the strong coupling limit, one has Ai ^ —l.bU and 
As ~ —3U from Sn — > 0.5, and A2 ~ —2U from the 
HF-RPA calculations[l3]. The corresponding coefficients 
are in good agreement with the DMRG results as can 
be seen from the inset. On the other hand, the results 
of HF-RPA deviate from the DMRG data apparently in 
the intermediate coupling regime, but this is quite un- 
derstandable since correlations involved in ([T]) cannot be 
accurately handled in HF-RPA when Ht and Hu become 
comparable, i.e. neither of them are perturbative. 

In summary, we have studied the low energy prop- 



erties of the one-dimensional half-filled 5/7(4) Hubbard 
model with the attractive on-site interaction by using the 
DMRG method as well as the perturbation theory. We 
found that the ground state is a CDW insulating state 
with the long range order in which the translational sym- 
metry is broken and all kinds of excitations are gapful for 
finite U < 0. Within our numerical accuracy, we found 
that the four-particle excitation gap is the same as the 
charge gap. While the charge gap (the four particle ex- 
citation gap) behaves non-monotonically, the others in- 
crease with increasing \ U\ and have a linear-U asymptotic 
behavior with different coefficients. Therefore, we believe 
that the one-dimensional attractive half-filled Hubbard 
model for the SU{4:) and SU{2) cases belong to differ- 
ent universality classes. Moreover, we find that the na- 
ture for the SU{4:) case can be further generalized to the 
other SU{N > 2) casesflSi]. At the end, it is worthwhile 
to mention that since the four-particle excitation gap as 
well as the charge gap are the smallest energy scale for 
the SU{4:) case with U <0, it would be very interesting 
to detect four-particle process (excitations) in an ultra- 
cold fermionic atom system with the hyperfine spin-3/2. 
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